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ABSTRACT 

With the growing interest in and ability of using weak tensing studies to probe the non-Gaussian properties of the matter density field, 
there is an increasing need for the study of suitable statistical measures, e.g. shear three-point statistics. In this paper we establish the 
relations between the three-point configuration space shear and convergence statistics, which are an important missing link between 
different weak lensing three-point statistics and provide an alternative way of relating observation and theory. The method we use also 
allows us to derive the relations between other two- and three-point correlation functions. We show the consistency of the relations 
obtained with already established results and demonstrate how they can be evaluated numerically. As a direct application, we use these 
relations to formulate the condition for E/B-mode decomposition of lensing three-point statistics, which is the basis for constructing 
new three-point statistics which allow for exact E/B-mode separation. Our work applies also to other two-dimensional polarization 
fields such as that of the Cosmic Microwave Background. 
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^1. Introduction 

The statistical study of the weak gravitational lensing effect by the large-scale structure, also called cosmic shear, is one of the 
established tools to investigate the matter distribution of the Universe and to constrain the parameters of cosmological models. The 
K*" undergoing boosting o f weak lensing survey size and quality gradually enables one to explore the small, non-linear scales where 
0\ _ a wealth of signal lies dBemardeau et al.L l2002t iPen et all [20031: iJarvis et al.L l2004t ISemboloni et al.L 1201 lb . Along with this trend, 
. more and more effort has been put into the development of statistical measures that can probe the non-Gaussian signal arising from 
the non-linear growth of stiTicture on those scales. Among them, the three-point (3-pt) statistics may be the most widely used. 

Two basic quantities considered in gravitational lens theory are the convergence k and the shear y. Defined as the dimensionless 
' surface mass density, /c is a weighted projection of the three-dimensional (3D) matter density contrast 6. The shear y, on the other 
' hand, is directly accessible from observations. Therefore, the theoretical framework of gravitational lensing should include the 
relation between configuration space k and y statistics as well as the one relating configuration space statistics to their Fourier space 
' counterpaits. At the level of two-point (2-pt) statistics, such relations have already been established. For 3-pt statistics, the relation 
, between the shear 3-pt coiTelation functions (}'3PCFs) and the co nvergence bispec t rum, w hich is the Fourier counterpait of the 

■ 3-pt convergence correlation function (/(•3PCF), has been derived bv lSchneider et al.l (l2005l) . The other non-trivial relation, the one 

■ between y3PCFs and /c3PCFs, is still missing. One purpose of this work is to establish this missing link. 

■ How to perform E/B-mode decomposition is also a major concern of the weak lei ising community . For observational data an 
E/B-mode decomposition provides a necessary check on the possible systematics (e. g. ICrittenden et aU, l2002t iPen et all |2002|) . In 
recent years there have been several efforts to construct better statistics which allow for an E/B-mode decomposition at the 2-pt 
level ([Schneider & Kilbinger, 2007; Eifler et al., 2010; Fu & Kilbinger, 2010; Schneider et al., 2010). They all use weight functions 
to filter the shear 2-pt coiTelation functions (y2PCFs), and the condition for E/B-mode decomposition transforms to a condition on 
the weight functions. Such a condition at the 3-pt level is also missing so far. We will see that with the aid of the relation between 
the 7'3PCFs and the /f3PCFs, one can easily formulate this condition. 

The paper is organized as follows: In Sect. 2 we show how the relation between the y3PCF and the /(■3PCF is obtained. In Sect. 3 
we investigate the coiTespondence between the derived relation and akeady established results. We then extend our results to other 
y3PCFs in Sect. 4, and in Sect. 5 we present an application of the 3-pt relations, deriving the condition for E/B-mode separation of 
3-pt shear statistics. How these relations can be numerically evaluated is demonstrated in Sect. 6, and we conclude in Sect. 7. 
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2. Relation between 3-pt y and k correlation functions 

2.1. The form of the relation 

At the 2-pt l evel, the relation between t he configuration space shear and convergence statistics is the ^+ 
dCrittenden et al.L l2002t [Schneider et all l2002h 



-I 



4 H(x -y) + 



^_ relation 



(1) 



where H and are Heaviside function and ID Dirac delta function, respectively. The functions ^+ and ^_ are defined as ^+(x) : = 
{kk}{\x\) and ^-(x) := {jy}{x) e"^*, with 0^ being the polar angle of the separation vector x, and ( ) indicating the ensemble 
average. Note that the shear y is a spin-2 quantity, i.e. it gets multiplied by a phase factor e"^"^' when the coordinate x rotates by 
(px- Consequently, (jjXx) has a spin of 4. Being the product of (yyXx) and a phase factor of e""*"^', the quantity ^-{x) no longer 
depends on the polar angle of x. 

The relation ([TJ has already taken both the statistical homogeneity and isotropy of the shear field into account and is therefore a 
one-dimensional relation of quantities on the real domain. The derivation of the ^+ -^^ relation originates from the relation between 
^+ and ^_ and the convergence power spectrum P^, 



PAt) = 27r I dxx Jo({x) = 2n I dx x 74(fo) . 
Jo Jo 

Inverting one of the relations in (|2) one can write ^+ and ^_ in terms of each other, e.g. 

r ^U{x)PJf)^ [ dyy^^iy) [ dl £ J^{{x) Mty) , 
Jo 27r Jo Jo 



(2) 



(3) 



and the final form of th e relation ([TJ can be reached by performing the ID Bessel integral whose result can be obtained from 
iGradshtevn eTaP (l2000l) . 

The same procedure, however, fails to work for 3-pt statistics since the corres ponding Bessel integra l actually consists of three 
integrals, and they have highly complicated dependencies on the arguments (see [Schneider et aU l2005l) . A brute force numerical 
evaluation of these integrals is also extremely challenging due to the oscillatory behaviour of the Bessel functions. 

Since the advantage of transforming to the Fourier plane and back no longer holds for 3-pt statistics, we attempt to stay in 
configuration space, which at least avoids the problem of oscillatory integrals. One can see from ([!) that the result of the Bessel 
integral in ([3]l is actually not oscillatory, as expected. 

The configuration space 3-pt shear correlator can be written as {y{X\)y{X2)y{Xi)), with Z,- being the positions on the two- 
dimensional (2D) plane where the shear signals are evaluated. Following the assumed statistical homogeneity of the shear field, 
the correlator depends only on the separations of these three positions. We choose xi = Xi - X3 and xz = X2 - X3 to be its 
arguments (see the leftmost sketch of Fig.[r|i and write the correlator as {yyy}(xi,X2). After the same procedure is applied to the 
3-pt convergence correlator, the relation we are interested in will be shown to be of the form 



{yyy}ixi,X2) = J ^^^^ J ^^^^ {KKK)(yi,y2) Go(xi -yi,X2 -yi) , 



(4) 



where we have defined the convolution kernel Go for which we need to find an explicit expression. 

Writi ng the relation in the form of a convolution is motivated by the Kaiser-Squires (K-S) relation between the convergence and 
the shear dKaiser & Squireslll993[) . 



dy K(y) D(x - y) , with the K-S kernel £)(z) = r 

2*2 



which yields the result ([U and also allows us to express the kernel Go as 

Goia,b)^-J dh'Div) D{v - a) D(v - b) ^ J 



d^v 



1 



1 



(v*-a*y (v*-by 



(5) 



(6) 



Here, for simplicity, we have adopted the complex notation for the K-S kernel, i.e. we have identified the 2D separation vectors 
with complex numbers. Throughout the text we will use the vector and complex notations interchangeably, and use x to indicate a 
complex quantity, x for its absolute value, and x* for its complex conjugate. 

The integral in ([6]l is difficult to perform directly, so we first take a look at the more studied 2-pt case. The relation between 2-pt 
y and k correlation functions can be written in the same way as 



(yyXx) 



d^y {KK}(y) F(x-y), 



with 



F(z) 



D(v) D(v - z) 



d\ 



1 

^2 



1 



(7) 



(8) 



(V* - Z*)2 ■ 

Unlike the case of the ^+ - ^_ relation, we have not assumed a statistically isotropic field for ([4]i or ([TJ. The ^+ - ^_ relation is 
actually what one should obtain after adding the assumption of isotropy to ([7). 
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2.2. The form of the convolution kernels 



Now we aim for obtaining the forms of the F and Go kernels, which can be seen as the 2- and 3-pt equivalence of the K-S kernel 

1 ^ "2 



Q. Introducing the symbols d = di + 182 and s 5? + c?? = dd*, the definitions of k and y read 



^=^VV, 7=^5V, (9) 

i.e. both the convergence k and the shear 7 are second-order derivatives of the deflection potential iff. It is then convenient to use iff 
as a link between k and y. Using the identities V In \x\ = and In \x\ = 2n6^^\x) which hold for a 2D x, one can easily verify 

the consistency of (|9]l with the relation between 1/' and k (e.g. lBarteknann & Schneiderll2001l) . 



,^(x)^^ J(fyK(y)ln\x-y\. (10) 

Applying the operator on both sides of ( fTOl l and taking (|9]l into account, one reaches the K-S relation (|5]), since !D(z) -d'^ In |z|. 
The same procedure can be generalized to second-order statistics. The 2-pt equivalence of ( fTOt is 

<^/^(xiMfe)> = ^ J d^yi \vi\xx-yx\ j d^y2 \n\x2-y2\ (^CyikCVi)) • (11) 

Using the statistical homogeneity of the k field, and re-defining the integration variables, (ITTt reduces to 

<(A(xi)i/r(x2)> = ^ 1 d^y {KK)(y) f d^u \n\u\ln\xi - X2 - y - u\ 

d^y {KK}(y) r'(xi-X2-y) , 



where we have defined 



J^d^M ln|H 



r'iz) = d^u ln|H|ln|z-M| . (13) 



Obviously, 'F' is infinite at every z, which is related to the fact that i/r is defined only up to an additive constant. However, we shall 
only need the derivatives of f^'. So we define 

riz) = r'{z)-r'm^ Jd^M in|H|inji^^J , (i4) 

and will use ;F and ;F' interchangeably. Let if denote the angle between u and z, (fT4l i can be rewritten as 

T (Z) - — \ duuinu \ di^ In 1 cosi^H — —\. (15) 

2 Jo Jo \ M I 

The integral over ^ yields zero if |z| < u, and 4;/rln(|z|/M) otherwise. Thus 

^(z) = 2;r P'dMM InM ln(|z|/M)= JlzP(lnlzl-l) . (16) 
Jo 2 

We are now ready to apply differential operators to (ITZt to get the relations of 2-pt shear and convergence statistics. As a 
consistency check, we first apply two operators to (fT2t . one acting on xi and the other on X2. According to (|9]l, this turns the 
l.h.s. of (fT2l i into 4 {k(xi)k{x2)}. On the rh.s. of (fT2l i the operators act exclusively on ;F, 

V2^^(jci -JC2 -y) = V2v2^(z) = V2(27rln|z|) = 4n^6'-^'>(z) , (17) 

with z = JCi - JC2 - 3' here. Using ( fTTl i. one easily sees that the rh.s. of (fT2l i after the operation gives 4 (a-zc) (jti - jci), which is 
equivalent to 4 {k{x{)k{x2)} under the assumption of statistical homogeneity of the k field. 

Now we apply the operator d\^d].^IA on (fT2l i. which turns the l.h.s. of ( fT2l l into (7(^:1)7(^1:2)). On the rh.s. the operation again 
acts only on T, 

}-dldlT(xx-X2-y) = ^d'Tiz) , (18) 



also with z = Xi - X2 - y. Remembering the definition of the kernel F O, this leads to 

■ z y 4 z 



F(z) = fd^ ^2 7T-^ = 7'^V(z) = 27r4T • (19) 
J v*'^ (v* - r 
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For the 3-pt kernel Go we split the integral in © into 
d V 



/ 

1 r 2 ± 

(a* - by j " v*2 



f 1 

[(V* -a*)2 

where we have assumed a + b. From (fT9l l as well as 



1 



(V - by 



{a* - b*)- 



f 



we obtain 



Go(a, b) 



I 



1 



— - — : 



1 



In 



(v*-a*)2 {v*-by {a* -by 



Z*2 

b_ 



V* - a* 



ft* 



2;r 



(a* - ft*)3 



b_ 



(20) 

(21) 

(22) 
b for 



The forms of the kernels (fT9l l and (l22t hold rigorously outside their singularities (at z = for F\ al a - Q,b - Q, and a = 
Go). One may wonder if additional delta functions exist at these singularities. We will show in Sect. |3]that this is not the case. 

The method we used to derive the forms of the kernels ([19) and (l22l i also allows one to derive the relations between other corre 
lation functions of weak lensing quantities in a systematic way. We present explicit forms of some of the relations in Appendix iBl 



2.3. The relations 

To summarize, we have obtained: 



and 



>2 {KKK}(yi,y2) 
yi -xi 



(y* -xy 

1 



(yi* - xi* - y2* + xzy \(yi* - xiy (yz* - X2*) 



yi-x2 



(23) 



(24) 



(y\* -xx* -^2* +X2y \Cyi* -x\y (j2* -J^i*)" 

In these relations we have applied the statistical homogeneity of the convergence field, but not the statistical isotropy. Making 
use of the latter, one can derive the ^+ - ^_ relation from (l2Jt . as will be shown in Sect. 



3. Consistency checks 

3. 7 . The case of uniform k 

There is a physical condition which will directly serve as a test of the k - y relations Q, ( l23T l and (l24l l. At the 1-pt level, for 
the K-S relation, a uniform convergence field does not result in any shear. At the 2- and 3-pt level, the physical condition could be 
that a uniform {kk) {{kkk)) field leads to a vanishing shear correlation (yy) ({yyy}). 

One can easily see that both the F and Go kernel we obtained satisfy this condition. If there are additional terms at the singular- 
ities of the kernels which contribute to the integral, a non-zero (yy) ({yyy}) term would be generated and the condition would not 
be satisfied anymore. Thus we argue that the expressions (l2Jt and (l24l i are already complete. 



3.2. Consistency witti ttie ^+ - f_ relation 

Now we consider whether (l23T l is consistent with the ^+ 
That the two relations are consistent is equivalent to 



- ^_ relation ([TJ, which can be regarded as the isotropic form of (l23T l. 



Jo 



d<p. 



y-x 



(y*-xy 

To verify that d25] t indeed holds, we attempt to solve the 

■ V — V 

d4>. 



4x2 _ I2y2 



H(x-y)- 



6(x-y) 



-integral on the l.h.s., 

-2n 



(y*-xy 



Jo 



d0 



y 



{y &-">' - xy 



(25) 



(26) 



where = 0,,- <Px has been defined. The 0-integral can be carried out using the residual theorem, yielding 27t{x^ - 3y^)/x'^ when 
X > y and zero when x < y. One can see that this result corresponds to the Heaviside function on the rh.s. of ( l25T l. 

At the singularity x = y the 0-integral is not well defined, which means one cannot rule out the existence of additional delta 
function at x = y in the result of the 0-integral. This ambiguity can again be eliminated by using the physical condition 'a uniform 
(kk) field leads to a null shear correlation {yy}', which translates to 'a constant ^+ yields vanishing ^_ ' here. In this case, a 
delta function is indeed required to satisfy this condition, and the prefactor of the delta function can be determined to be n/2x, in 
consistency with ( l25T l. 
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3.2.1. K-S relation and its isotropic form 

A similar consistency exists between the K-S relation and its isotropic form. As both forms are already well-known, they can 
serve as a further support for our argument. 

For an axisymme tric distribution of ma tter, i.e. k(x) = k(x), the following relation is established between the shear and the 
convergence (see e.g. lSchneider et al.lll992h 

y(x)^[K(x)-K(x)] e^'^' , (27) 

with K defined as 

^(^):=4 r ydyKiy). (28) 

X Jo 

This is equivalent to 

y(x) = - ^ jydy K(y) [2H(x -y)-x6'^\x- y)] . (29) 

In the case of a uniform convergence field k{x) - const., one can see that the integral of the Heaviside function and the delta function 
parts cancel each other 

The similarity between ( |29T l and the ^+ - ^_ relation ([TJ is remarkable: they both have integrals of a Heaviside function part and 
a delta function part which cancel each other for constant k and {kk), respectively, and the 2D correspondences of both do not have 
an additional delta function at their singularities. 

3.3. Fourier transformations 

In the Fourier plane the relation between the shear and the convergence has a simple form. With g denoting the Fourier trans- 
formation of the quantity g, it reads 

7(0= e^^Kit), foxti^Q, (30) 

with j6 being the polar angle of the wave vector t. The identity ( l30b can be obtained by relating the Fourier transforms of the 
definitions (|9]l of y and k. It directly reflects the fact that y is spin-2 while k is spin-0, and leads to the well-known result Py - P^. 
Comparing dSOl l to the configuration space relation (|5) which means y = /c * D/;r, using the convolution theorem, one can see that 
e^^ = D/n for ( ^ 0. 

The Fourier plane correspondences of (l27t and (l24l are also readily obtainable from the identity (l30l l. as 

<77)(0= t'^{rK)(f), forZ-^O, (31) 

and 

{yyy)ih, h, h) = s^'^p^^p^-^p^-^ {nmeu h, h) , for d, (3*0, (32) 

with /?, denoting the polar angle of d. These equations show that e^'/^ = F(e)/^^ for t * 0, and that t'^^^^+h+Pi) = -Go(^i, ^2)/^^ 
for (i ^ (2 and {3 - -ti - ti +Q, since Fjn^ and -Go/n^ are the convolution kernels for the configuration space relations by 
their definitions. 

In AppendixlAl we show explicitly that the Fourier transforms of F/tP' and -Gq/tt^, with F and Go given in ( fT9] ) and (l22l) . are 
indeed e*'^ and e^^^'^P-'^P^^ respectively. However one cannot obtain the forms of F and Go kernels simply through inverse Fourier 
transforming the phase factors e'^'^ and e^'*'^'^^-^^^*. This is due to the fact that the Fourier inversion theorem is valid strictly only 
for square-integrable functions, which is not the case for the phase factors. The same situation occurs for the K-S kernel D. 



4. The other ySPCFs 

Until now we have considered only the 3PCF of shear itself (7(Zi)7(Z2)7(X3)), which is one of the four independent pos- 
sible combinations considering that 7 is a complex quantity. The other three are (7*(Zi)7(Z2)7(X3)), (7(Zi)7*(Z2)7(X3)), and 
(7(Zi)7(Z2)7*(X3)), according to the choice made in fSchneider & Lombardil (l2003l) . Following [Schneider etaLldBol . we denote 
these four 73PCFs by r^'^j,(jci, JC2) (/ = 0, 1,2,3), with 'cart' emphasizing that the shear is measured in Cartesian coordinates, 
^cart = (TTT)' and r^'^j, (; - 1, 2, 3) corresponding to the 73PCF with 7* at position Z,-. Since we have considered statistical homo- 
geneity of the shear field, the Fcait's depend only on the separation vectors of the position Zi, Z2, and Z3. The other 73PCFs, i.e. 
those with two or three 7*'s, can be obtained by taking the complex conjugate of the Fcait's. 

Note that r'-JJ^{xi, X2) = <7*rr>(-Ki,-K2), rf^,^{xi,X2) = <7r*r)(-Ki>-»^2), and r'^^^^(xi,X2) = <7rr*>(-»^i.-»^2) are different functions, 
since Xi (X2) is defined to be the difference of the positions of the first (second) and the third 7 in the bracket. Due to the same reason, 
they can be transformed into each other through permutations and flips of the vertices of the triangle formed by their arguments (see 
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Fig.[TJ, and thus are not independent if argument permutations and flips are allowed. As an example, one has 





- ^ cart'--*-l'-*-2^ 




/~.*^v \~.^v \~.^v \\ 






= <r(Z2)r (Xi)r(X3)> 






= (r(Z3)r*(Xi)r(Z2)) 






= <y(Z2)r(Z3)r (Zi)) 


= rfartfe - Xl, 


-^1) 


= <r(Z3)r(Z2)r*(^i)> 


^I<^J^(-XuX2- 


-Xl), 



(33) 



where different lines correspond to different ways of labeling the same triangle with side lengths xi, X2, and \xi - Xzl- The same 
permutations and flips also reveal the inherent symmetry of rfj^, 

rS*„(-^l, ^2) = rfjjx, - X2, -X2) = Tfl{X2, Xl) = Tfl{-X2, Xl - Xz) = Y^Z{X2 - Xl, -Xl) = lt!j-Xi,X2 - Xl) . (34) 




In the case that the shear is measured re l ative to a center of the triangle, F^'^^j transforms to P'\ the natural components of the 
73PCF as defined in lSchneider & Lombardil(l2003l) . For a general triangle configuration, all four Fcan's are expected to be non-zero, 
thus all of them should be used to exploit the full 3-pt information of cosmic shear. 

Before relating the other Fcan's to the a-3PCFs, we extend /c to a complex quantity k = + ik^. Although the physical con- 
vergence is a real quantity, the convergence field corresponding to the measured shear signals can have an imaginary part due to 
e.g. systematical errors and noise. The shear component which corresponds to this unphysical imaginary part of the convergence 
field is identified as the B-mode, on which we will elaborate more in Sect. |6] When taking the B-mode into consideration, the 3-pt 
correlation functions of the convergence field can be written as 



(35) 



Apart from the E-mode (k^k^k^) term, there are still additional B-mode contributions to the real parts of the K's, namely (k^k^k^), 
{k^k^k^}, and {k^k^k^}. T he imaginar y part of the K's are composed of the parity violating terms which are expected to vanish due 
to parity symmetry (ISchneideri |2003|) . The property of K^''' under permutations and flips of the vertices of the triangle formed by 
their arguments is the same as that of F*'^j^. 

Similar to dU, the relations between F^^ and K^'^ for / = 1,2,3 can be written as 

r';}^(xi,X2)^~ JdV JdV^*'*(3'i,j2)Gi(;t:i-3-i,X2-3'2), (36) 



and 



rf^(xi,X2)^~ Jd^yi j'd^y2K'^^^(yi,y2)G2(xi-yi,X2-y2): 
^fJa(^i,X2) f d'yi J K''''(yi,y2) G^ixi -yuX2 -3-2) . 



(37) 
(38) 
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where the convolution kernels Gi, G2, and G3 have been defined. Again with the aid of the K-S relation, we can write these 
convolution kernels as 

Gi(a, ft) = - r d^v D{v) &{v - a) D{v - b) ^ fd'v ^ , , \.,2 

J J v*2 (v - ay (v* - by 

J ^ v2 (v* + a*)2 (v* +a* - by ' 

G,ia, ft) = - 1 d^v !D(v) Div - «) D*(v - ft) = J d^ ^ 

= fd^ 1 

J V- (v* -a* + by (v* + by 



(39) 



(40) 



and 



G3(a,ft) = - r d^v £)*(v) £)(v - a) £)(v - ft) = f d\' ^ — 

J J (v* - a* 



y (V - by ' 

When a b, the product of the three terms in the integrand of (|411 can be split into products of two, as 



1 1 



1 1 



1 1 
+ 



(V* - ay (V* - by 



1 



(a* - by v2 



V* - a* V* - ft* 



(41) 



(42) 



v2 (V* - ay (V - by (a* - by v2 
These terms are also obtainable from doing derivatives to the kernel T', 

J (v* - ay 4 

r d^vi^^ V'w=-. 

J V* - a* 2 a 

This way we obtain the form of the convolution kernel G3. The forms for the kernel G\ and G2 can be obtained likewise. The results 
are 



(43) 



(44) 



c,(».») = ^K'c).«w-.)]-|^(i.^), 



(45) 
(46) 
(47) 



The symmetries in the Fcart's and K's are also reflected in the G kernels. One can verify that G2(a, ft) = G\{b-a, -a), Gi{a, ft) = 
G\{-b,a-b) as results of the symmetry under permutations, and G2(a, ft) - G\{b, a), G^ia, ft) = G3(ft, a) as results of the symmetry 
under flips, in consistency with (l33T l. Similarly, one has Go(a, ft) = Go(ft - a, -a) - G^i-b, a - b) - Go{b, a), in consistency with 

01. 



5. Inverse relations 

So far we have obtained the expressions of the four y3PCFs as functions of the /cSPCFs. Written in a short form, they are 



1 



where ; runs from to 3. The forms of the G, kernels are given by (l22l i. (l45T l. (l46T l. and (l47t . 
These relations can be inverted. We define the kernels of the inverse relations to be G', i.e. 



<0 



,3 / cart ■ 



Using the convolution theorem, it is apparent from (|48) and (|49l > that 



3 ' 



From the corresponding Fourier plane relations of d48l . we also know 



Go 

7r3 



" ^3 



G2 
" ^3 



,2iOS,-ft+ft) 



G3 
" ^3 



,2it8,+/32-A) 



(48) 



(49) 



(50) 



(51) 
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which implies 

GoGq — GiG\ — G2&2 — G3G3 = 71^ . (52) 

Comparing (ISOl l and (|52] | one has 

g; = g; , (53) 

and further, 

g; = g; , (54) 

i.e. the convolution kernel for the inverse relation is the complex conjugate of the original kernel. 

This property of the convolution kernel has its root in the fact that y and k differ only by a phase factor This fact also endows 
the convolution kernels in the 1-pt and 2-pt relations between y and k with the same property. As is well known for the 1-pt relation, 
the inverse relation of (2) is 

K{x)^^ ^ A^yy(y)D*(x-y), (55) 
where the kernel is the complex conjugate of the K-S kernel D. The inverse relation of the 2-pt relation d?) can also be shown to be 

{kk){x) ^^Jdb {77}(y) F*(x - y) . (56) 



6. Condition of 3pt E/B decomposition 

Being mathematically a polarization field, the shear field can be decomposed into a curl-free component and a divergence- free 
component, usually called the E-mode and the B-mode, respectively. Performing such a decomposition when treat ing cosmic shear 
data has long been rec ognized as a necessity, since it provides a valuable check on the possible systematics (e.g. ICrittenden et all 
l200llPen etani200l . 

The E/B-mode decomposition can be done either on the shear field itself (e.g. lBunn et al 1 120031: iBunn'. '201 Oh. or at the level of 
correlation function s (e.g. Schneider, 2006). The complex survey geometry after masking, which is especially characteristic for a 
lensing survey (e.g. lErben et all |2009|) . renders the first option barely feasible, and singles out the correlation function as the basic 
statistic to be applied directly to the data. Thus the natural way to perform the E/B-mode decomposition on cosmic shear data is to 
derive statistics based on the shear correlation functions. 

A commonly used statistic for this purpose is the aperture mass statistic which can be expressed as a linear combination of ^+ 
and ^_ (.Schneider et al.>.200Z) . 



1 r 



~^ 



^+(i?)r7l-)-^_()?)r!P 



(57) 



where the forms of the weight functions T^^ and Tt^ are given explicitly in lSchneider et al.l(l2002b . The chosen forms of the weight 
functions guarantee that (^^ap) responds only to the E-mode and ( m^^ only to t he B-mode. 

The aperture mass statistics has been generalized to 3-pt level bv lJarvis et al ] (12004 and lKilbinger & Schneided (l2005h. and is the 
only st atistics available up to now which allows an E/B-mode decomposition at the 3-pt level. However, as found bv lKilbinger et al.l 
(1200 6'). it cannot ensure a clean E/B-mode decomposition when applied to real data. The lack of shear-correlation measurements 
on small and large scales, which arises from the inability of shape measurement for close projected galaxy pairs and the finite field 
size, prohibits one from performing the integral in (l57i from zero to infinity, and thus introduces a mixing of the E- and B-modes. 

In recent years, there have been several eff'orts to construct better statistics which allow E/B-mode decomposition 
dSchneider & Kilbingeit 120071: lEifler et al.L 120101: iFu & KilbingeR 1201 Ot [Schneider et all 1201 (A . all of them focusing on the cos- 
mic shear 2-pt statistics. These new statistics are based on the idea that the weight functions T^^ and Tl^ used in the aperture 
mass statistics are just one ex ample out of the many possibilities. In general one can define second-order statistics in the form 
dSchneider & Kilbingeit l2007l) 

Xoo 
d§ [^+(j?)r+(i?) + ^_(i?)r_(j?)] , 

^00 (58) 

BB = I d-d§ [^+id-)T+(&) - ^_(i?)r_()?)] , 
Jo 

for which the condition that EE responds only to E-mode and BB only to B-mode is found to be 



I § d§ T4^)Jo{m) = I )?dj?r_()?)74(^t?), orequivalently 
Jo Jo 



" 4 12.=- 



T+(§) = r_(i?) + 9 do T^e) -r 



02 04 
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Note that instead of being functions of the separation length as the aperture mass statistics, EE and BB are just numbers. At first 
sight one seems to have reduced the information quantity by integrating over the scale dependence in (ISSl l. In fact, the information 
can be easily regained by constructing a set of weight functions satisfying ( |59] |. As one example, (^^p) (0) and (a^I) (0) for any 
value can be reconstructed in the framework of (|58) by specifying T^(§) = {§/e) 1 6^ and r_(j?) = r!'' (§16) 19^. 

Since the condition for E/B-mode decomposition ( |59] ) still leaves large freedom for the choice of the weight functions, one can 
construct statistics which fulfill additional constraints, e.g. a finite support over the separation length. If one requires that EE and 
BB respond only to ^+(??) and with i^min < '& < )?max, where j^min and j^^^ax are the chosen small- and large-scale cutoff, T+{&) 
and must vanish outside the same range. Since T+ and are interrelated by (|59^ . one can specify only one of them to satisfy 
this constraint. The requirement that the other weig ht func t ion als o vanishes outside the specified range needs to be put as additional 
integral constraints. As shown bv ISchneider & Kilbingerl (|2007h . if one chooses T to vanish for {)■ < i^min and > J^max, then to 
allow an E/B-mode decomposition on a finite interval ??n,in < & < §max, T has to satisfy additionally. 



1. T"--''^-"-]. 51 



(60) 



which would guarantee that T+ vanishes for )? < -d-jnin and > i?max- 

Similar statistics are needed at the 3-pt level as well. The first step required is to formulate the conditions for 3-pt weight 
functions to allow E/B-mode decomposition, in analogy to ( |59] |. As we will show in this section, the relations between the y3PCFs 
and a:3PCFs that we derived provide a natural way of formulating such conditions. 

A pure E-mode shear 3-pt statistics is related only to the E-mode /(■3PCF (^k^k^k^^ but not to other 3PCFs with contribution. 

Therefore we first write the 3PCFs of and as linear combinations of the real and imaginary parts of the /r*'*'s, using i35[ . and 
then relate them with the Fcart's through (|49T l, as the Fcan's are the directly measurable statistics from a lensing survey. The results 
read 



{k^k^i^} = -Re 



1 

-Re 
4 



(K^i^K^) = ^Re 



and 



/ B B B\ 
{k K K ) 



1 

-Re 
4 



Im 
Im 
Im 
Im 



_/^(0)_^(l)+^(2) 



^Re [ G; * F<°', + Gl * T^^ + * Ff,', + G; * ] , 

[-g; . rZ - Gl * rZ + * i<Z + g; * f« ] , 
[-g; * rZ + Gl * rZ - Gl * F^', + g; * f« ] , 
_^(0) ^ ^(1) ^ ^(2) _ ^(3)j ^ , rfj, -H Gl * rli> -H g; * f^^ - g; * f« ] , 



_/^(0) + /^a)_/^(2) + 7r(3)j 



(61) 



7r(0)_7r(l) + /^(2) + ^(3)j 
7^(0) +7^(1)- 7^(2) + ^(3)] 
7^(0) +7^(1) + 7^(2). ^(3)] 



4;r3 
1 

1 

'4^^ 



GS 


* i cart 


-Gt 


* ^ cart 


+ G5 


*F<2) 

cart 


-hg; 


, r(3) 

''^cart 


GS 


*F<"> 

cart 


+ GI 


*F<'' 

cart 


-G5 


*F<2' 

cart 


-hg; 


cart 


GS 


*F<"> 

cart 


+ GI 


*F<'' 

cart 


+ G2 


*F<2' 

cart 


-g; 


*F<3' 

cart 



] 



1 



(62) 



which shows how the E- and B-mode /(■3PCFs can be computed when the full information of the Fcart's is available. In the ideal 
case that there exists no noise or systematical effects, only the E-mode term k3VC¥ ^/c^/f^/c^^ is expected to be non-zero, since it 
corresponds to the correla tion in the physical density fiel d which leads to the correlation in the shear signal. 
Following the ideas of lSchneider & Kilbingeii (l2007h . we construct a new statistic 



EEE= J d^xi ^ (fx2 {k^i^j^) {xi,X2) U{xi,X2), 
which by definition responds only to E-mode. With the help of ( I6TI 1 we can link EEE to the observable Fcart's, as 



(63) 



EEE 



d X2 U(xi,X2) Re 



1 

4^ 
1 

4^ 



Re 



Re 



I dVi I d^yi G*ixi - yi,X2 - y2) r'il^iyuyi) 
J dV J jF«„Cyi,3-2) J d^xi J (fx2 U(xuX2)G.(xi-yi,X2-y2 

j dVi j d^yi 2rltCyi>J2)(G**f/)(ji,j2 



(64) 
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where in the first equation we have specified f/ to be a real function, and in the second equation we have used the fact that 
Gi(-a, -b) - Gi(a, b). 
Denoting 

G**U-. r« , (65) 

the expression of EEE (|64] | has a similar form as (ISFt . We can see in this form that BEE responds only to the E-mode if the weight 
function T's satisfy 

* Go = t/ = * Gi = T^^^ * G2 = T*-^^ * G3 . (66) 

with U being a real function. One can easily verify that these conditions are also necessary conditions. Noticing that T^'^ * G, is the 
corresponding weight on K^'\ the condition that these functions being purely real is required to separate the parity-violating and 
non-violating terms in (|35] |. In addition, ( |66] | is required to cancel the parity non-violating B-mode terms {k^k^k^}, {k^k^k^} and 

The statistics containing the contribution from only one of the B-mode terms can be constructed in the same way. Omitting the 
arguments for notational simplicity, they can be expressed as 



) r('' 

cart 



EBB =-Re\ ff Y^T^' 

BEB=i/?J fT^^^r^ 

L^^ !=0 



BEE 



and 



BEE =-/m 
4 



EBE 



EEB 



,=() 



) r('' 

cart 



) r('' 

cart 



;■={) 
3 



BBB 



) r('' 

cart 



with r^°> * Go = * Gi = -r<2' * G2 = -t^^^ * G3 , 
with * Go = -r*" * Gi = r<2' * G2 = -r<^> * G3 , 
with * Go = * Gi = -r*^> * G2 = * G3 , 

with * Go = -r^" * Gi = r^^' * G2 = r^^' * G3 , 
with * Go = r<'^ * Gi = -r^^^ * G2 = r^^^ * G3 , 
with r*"' * Go = * Gi = r*^* * G2 = -r^^' * G3 , 
with r*"' * Go = -r<" * Gi = -r(2> * G2 = * G3 . 



(67) 



(68) 



For all of them the condition that T*"' * Go is real has been imposed. 

The four parity violating statistics (l68t can be used as a check on parity violating systematical errors, while the other B-mode 
statistics ( l67b allows for a further examination of the B -modes. 

With ( l52b one can easily invert the conditions on the weight functions to express the weight function T's directly in terms of 
each other Take the conditions for EEE (l66t for example. One can write the Fourier transforms of T^^\ T^^' and T*^* as functions of 
the Fourier transform of T*"* as 



f (1) = J_f (0) GqG* fP) = IfmcoG;, f = If (0) Gog; . 

7T^ jfi 



(69) 



In order to simplify these relations, we now attempt to give simple expressions for GoG* for / = 1,2,3. With (ISTT l one has 
Go{{i,t2)G\{ti,t2) = e"*'^', which does not depend on {2. Noticing that F = tt^ e'*''' (see Sect.[33|, we actually have 
Gditi, t2)G*Ati,t2) - n'^F(ti), which yields in configuration space 



iGo*Gl)(a,b)^K'F(a)S'^\b). 

Similarly one can derive that 

(Go * GI) (a, b) = n^Fib) 6^^\a) , (Go * G;) (a, b) = n''F{a) 6^^\b - a) . 
Inserting dTOl i and ( fTTT i into (|69] l, one obtains 

T^'\xi,X2) \^^y T'^'^iXl -y,X2)F{y) = ^ J t'^^)(xi - y,X2)j-^ , 

T^'-\xi,X2)^\ r d^yT^''\xi,X2-y)F{y)^- f t'-^\xi,x2 - y)^ , 
J n J y'^ 



(70) 
(71) 

(72) 
(73) 
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T'^Xxi^Xi) = J d^y -y,X2- y)F(y) = ^ J d^y -y,X2-y)j;^, (74) 

where we have used the expression of F ( fT9l l. 

To summarize, if one chooses an arbitrary form of r*°^ which makes r'"* * Go real, constructs T^'*, T*^' and T*^^ according to 
(l72l l. ( l73l l and (l74l l. and uses these weight functions to weight the measured r[,'^'s, then the resulting statistic EEE as defined in ( l64l i 
receives contributions only from {k^k^k^) but not the terms affected by the B-mode. The B-mode statistics can be obtained from the 
measured T'^I^'s through a similar procedure. Equations (l72t-(l74l are the analogue of (l59b for 3-pt functions. 

We note again that T''^^^^, T*^^^ and F^'^j are not independent under transformation of their arguments. Thus the statistics EEE ( l64l i 
as well as the B-mode statistics ( |67] ) and ( |68] ) can all be written in terms of linear combinations of r*°'^ and F''^^ alone. However we 
shall keep the current redundancy since it allows for simple analytical expressions of the relations between the weight functions. 

So far one still has much freedom in choosing the form of T^^^\xi,X2). This freedom can be exploited to construct statistics 
which do not respond to the y3PCF at smaller or larger angular separations than can be probed by the survey. We leave this to future 
work. 



7. Numerical evaluation 

7. 1 . Design of the sampling grid 

We have written configuration space relations between weak lensing statistics in the form of convolutions, e.g. (l23T l and ( l24l i. 
where the convolution kernels are complex, have non-trivial spin numbers, and feature singularities. The convolutions can be per- 
formed numerically, but special care must be taken of these properties of the integration kernel. 

One can take the K-S kernel as an example of this kind of convolution kernel. The K-S kernel has an integer spin of 

2, so an azimuthal integration of the kernel around its singularity at z = should give zero, i.e. the values of the kernel along the 
circle cancel themselves due to their opposite phases. This property renders its singularity harmless, but entails the condition that 
the sampling grid should guarantee the cancellation. Such a requirement of a special grid design has already been realized in early 
lensing mass reconstruction works (e.g. ISeitz & Schneiderl [l996). For a spin-2 kernel like the K-S kernel, a common square grid 
already suffices if the singularity is placed at a center of rotational symmetry, i.e. either onto a grid point or at the center of four grid 
points. In the case of the former, the grid point at the singularity has to be discarded. The latter, as shown by the left panel in Fig.|2] 
is a better choice considering the sampling homogeneity. 



spin 2, 3 



spin 2, 3, 4 



spin 2, 3, 4 



Fig. 2. Examples of sampling grids applicable for 2D integration over singular kernels of different spin values. The cross in the 
center indicates the position of the singularity in the integration kernel. 



In general we need to deal with convolution kernels with different spin numbers. For example, the kernel F between y2PCF 
and /c2PCF ( fT9] l. which is proportional to z/z*^, has a spin of 4. In this case the square grid cannot guarantee the phase cancellation 
around the singularity any more (see Fig. |3). With a similar analysis as shown in Fig. [3] one can see that a triangular grid (middle 
panel, Fig.|2ll can actually guarantee the phase cancellation around the singularity of a spin-4 kernel. When using a triangular giid, 
the singularity can also be put either on a grid point or at the center of three giid points. The former choice loses the grid point at 
the singularity, but is applicable to spin-3 kernels where the latter fails. 

To achieve a high numerical accuracy, it is required that the circle integral around the singularity is well-sampled. If one uses 
a square (triangular) grid, the innermost circle is only sampled by four (six) grid points, which is not enough for many /(--models. 
To remedy this, one can duplicate the sampling grid, rotate it around the singularity, and put it on top of the original grid. We show 
the result with the square grid and 45 degrees of rotation in the right panel of Fig.|2] The resulting grid is also applicable for spin-4 
kernels, as shown in the right panel of Fig. |3] This non-standard way of constructing sampling grids, although not creating the best 
grids in terms of sampling efficiency, deals very well with the singularity of the integration kernel, and can easily generate sampling 
grids applicable for kernels of any spin number We use this kind of grid in our numerical sampling. 
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spin 2 spin 4 



Fig. 3. A square grid guarantees phase can- 
cellation around the singularity of spin-2 ker- 
nels (left panel), but not that of spin-4 kernels 
(right panel). The crosses indicate the posi- 
tion of the singularity of the integration ker- 
nel, while the dots are the grid points closest 
to the singularity. The polar angles of the ker- 
nel at the grid points are indicated by the di- 
rections of the arrows. For a spin-2 kernel they 
cancel each other on a square grid already. For 
a spin-4 kernel they cancel each other if the 
square grid is duplicated, rotated 45 degrees 
and put on top of the original grid. 



Additional complications arise when performing the integration for 3-pt statistics. There are two 2D integrals in (@), ( l36l l. 
(l37l i and ( l38l l. The corresponding integration kernels (l22t . ( |45] ). ( |46] | and ( l47l l all have three singularities. Luckily we can split 
each integration kernel into four additive terms and perform the integrals over each of them separately. Moreover, one can apply 
translational shifts to the integrands so that in each 2D integral there is only one singularity in the integration kernel. The singularity 
can also be shifted to the origin of the grids for numerical simplicity. After all these procedures, the four relations can be written as 



riart(-«^i'-«^2) =^ Jd^yi J d^y2^-[{KKK}(y2+xi,yi+y2+X2) + {KKK)(yi+xi+y2,y2+X2)] ^ 

1 V2 1 

+ [{KKK)(y2 + Xi,yi +y2+X2) + {KKK){yi+Xi +y2,y2 + X2)] — r — r \ , 

3'! yi^j 

r*2t(-ti,-C2) ^ ~ ~ J '^^y i{K*KK}{xi,y + X2) + {K*KK}(y + xi,y + X2)) ^ 

2 C C 11 
-H — d^yi d^y2 ({K*KK}iy2 + Xi,yi + X2) - (k* KK)(yi + y2 + Xi,yi + X2)) —r— , 
J J yi ^yi 

\ r 1 

^cL^xi,X2) = - - \ d^y {{KKK)(y + xi,y + X2) + {KKK)(y + Xi,X2)) — 
Tll^{x-s,,X2) = ~ ~ j ^('^'^'^ )(xi,y + X2) + (kkk }(y + xi,X2)) 



2 1 1 

d y2 i{KK*K)(yi + xi,yi +y2 + X2) - {KKK)(yi + xi,y2 + X2))—t— , 

yi*^ y2 



7 1 1 

d y2 ({kkk'Xxi +y2,yi+y2 + X2) + {KKK){yi +y2+Xi,X2 +y2)) —;— ■ 

yi*^ yi 



(75) 



(76) 



(77) 



(78) 



We can see that the integration kernels are either spin-2, spin-4, or spin-3. All sampling grids shown in Fig.|2]are applicable to spin-3 
kernels. 



7.2. Numerical results for 2PCF 

We now construct several toy models for 2-pt and 3-pt convergence and shear correlations, and use them to test the relations 
derived as well as the numerical sampling method. 

In the 2-pt case, we build two models for the convergence correlation function: (kk) (r) - 1/ r, and {kk} (r) = e"*^ . Using the 
well-established ^+ - ^_ relation ([T) we can obtain the corresponding models for the shear correlation function: (yy) (r) - e"*'*' /r, 
and (yy) (r) = e'*''^' [(r^ + + 6) + 2r^ - g] /r^. 

Fig.|4]shows the comparison between these shear correlation function models and the shear correlation functions sampled using 
( I23] ), with the corresponding convergence correlation function models as input. The numerically sampled values match the analytical 
models very well. 



7.3. Numerical results for 3PCF 

We build models for 3pt shear and convergence correlation functions via the 3pt correlation function of the deflection potential 
t//. Suppose we evaluate the fields at positions X, Y and Z. By definition we have 

{KiX)KiY)KiZ)) = (j^i) (5 Vf) (^Vij <^(Z)«A(10(^(2)> , (79) 
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Fig. 4. Left panels: {yy){r) with fixed 
(pr = 7r/6 as a function of r. Right pan- 
els: Polar angle of (jyXr) with fixed 
r = 20 as a function of (pr- The black 
curves are the expected values while 
the dots are the results from numeri- 
cal evaluation of 2D integral in (l23l l. 
The (kk) used in the upper panels is 
{KK){r) - l/r and in the lower panels 
{KK){r) - exp(-r-). 



V 









Re, analytical 

Re, numerical 

_ Im, analytical 

» Im, numerical 
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Fig. 5. Comparison of numeri- 
cally evaluated (yy/Xx,)!) and 
{yyy*}(x,y) with their analytical 
toy models described in this sec- 
tion adopting a = 0.05. 'Re' and 
Tm' indicate the real and imag- 
inary parts of the shear correla- 
tion functions. Left panels: The 
functions are evaluated at x - 
0.25k e-'"'^ andy = O.llk e'"'^ 
for different values of k; Right 
panels: The functions are evalu- 
ated at jc = 0.75 e^'"'^ and y - 
0.45 e"^ for 50 equally spaced cp 
values. 



and 



{y(X)y(Y)y(Z)} = (^^ll^^ll^^zl {<p(XMY)ifr(Z)) , 



Now we assume a model for {iJ/(X)i//(Y)il/(Z)) as 

<^(Z)iA(IO<A(Z)> 



1 

8c? 



(80) 



(81) 



(82) 



with X - X - Z and y - Y - Z. This model is special in the sense that it does not depend on the angle between x and y, but is 
nevertheless simple and rather local. The statistical homogeneity of the field enables us to write the 3pt correlation function as a 
function of two 2D spatial coordinates, which we have chosen here to be x and y. 
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Performing the derivatives, we find the corresponding 3pt shear and convergence correlation functions also depend only on 
X - Z - X and Y - Z - y, and read 

-Q-(:t^+/) 

{KKK){x,y) = [ - 4 + Of {(i{x^ + /) + 8jc • j) - a^i^x^ + y^f + 4(x^ + y^)x ■ y + (,x^y'^) + a'x^y^\x +yf\ , (83) 

{jjj){x,y) = x^y^ix + yf e-^'-^''^-^'^ , (84) 

and 

{yyy*){x,y) = [ix^ + %xy + 2y^ - 4-axy \x + y\^ + a^x^y^(x* + y'f] e-"**'+'''' . (85) 

Using ( l83T l as the input model for (kkk), we numerically evaluate {'yyy)ix,y) and {'yy')'*){x,y) using (l75l l and ( iTST l. The results 
are then compared to the analytical models for the 3pt shear correlation functions (l84l i and (ISST l. As shown in Fig. |5] the numerical 
evaluations closely match the analytical models. 

Hence, we have proven numerically that the relations between y3PCFs and a-3PCFs ( fTSl ) and ( iTST i are correct, no additional 
delta functions are needed. At the same time we have shown that these relations can be numerically evaluated with a high accuracy. 
Therefore these relations provide a better way of relating the observable shear signal to the underlying matter density field than the 
original way, i.e. using the relations between the 'y3PCFs and the k bispectrum, since the latter does not allow for an easy accurate 
numerical evaluation. 

8. Conclusion 

We derived the relations between the 3-pt shear and convergence correlation functions which had been an important missing 
Unk between weak lensing three-point statistics. As an intermediate step, we found the 2-pt analogue of these relations and proved 
that it is the non-symmetrized form of the existing ^+ - ^_ relation. By drawing analogy to the corresponding 1-pt relations, namely 
the Kaiser-Squires relation and its isotropic form, we have further revealed that the newly derived relations and already established 
results fit into the same theoretical framework. The consistency of the configuration space relations with the known Fourier space 
relations have also been shown. 

The 3-pt relations derived are simple both analytically and numerically. They can be used as an alternative way of relating 
the measurable 3-pt weak lensing statistics with the statistics of the underlying matter density field. Up to now one has to use 
the relations between the y3PCFs and the convergence bispectrum to link theory to t he observable 3-pt she ar signal. Since the 
73PCFs are very oscillatory and complicated functions of the convergence bispectrum (ISchneider et al.Ll2005l) . it is hard to study 
the behavior of the 3-pt shear signal for a given convergence bispectrum model. With the relations we derived, one can instead study 
the properties of the 3-pt shear signal by constructing models for the /c3PCFs. 

The method we used to derive these relations is based on the relation between the 2-pt correlation functions of the lensing de- 
flection potential and the convergence. The same method also allows one to systematically derive the relations between correlations 
functions of other weak lensing quantities, including the deflection potential t//, the shear y, the convergence k, and the deflection 
angle a. We present the forms of some 2- and 3-pt relations in Appendix. IbI Some of them are potentially of interest for studies of 
galaxy-galaxy(-galaxy) lensing and lensing of the Cosmic Microwave Background. 

Since the relations we obtained have complex kernels with non-trivial spin number and singularities, special care is needed when 
they are used numerically. We demonstrated how the numerical evaluation can be done, in particular the design of the sampling 
grid. Examples of numerical evaluation were shown for both 2- and 3-pt relations using toy models for the convergence correlation 
function. Their results match very well with the analytical expectations. 

Separating E- and B-modes from measurements of the y3PCFs is particularly important since the systematic effects at the 
3-pt level are less understood. So far the only 3-pt statistics allowing for an E/B-mode decomposition is the aperture mass statis- 
tics dJarvis et al.', '2004^. 'Schneider et all l2005l) which is plagued with the same problem as the 2-pt aperture statistics pointed out 
by ^ilbinger et al. (2006). To amend t his problem, one needs to construct the 3-pt correspondence of the newly developed 2-pt 
statistics (Schneider & Kilbinger. l2007l: lEifler et al.L |20 1 Ot IFu & K ilbingei', '20 1 0!: ISchneider et all |20 1 0) which allows for an E/B- 
mode decomposition on a finite interval. As a direct theoretical application of the 3-pt relations derived in this paper, we used 
them to formulate the conditi ons for E/B-mode decomposition of lensing 3-pt statistics, in analogy to the 2-pt condition given by 
ISchneider & Kilbingeii ( l2007 l). These conditions are the basis of formulating additional constraints which lead to E/B-mode decom- 
position over a finite region, therefore they provide a starting point for future works on constructing better 3-pt shear statistics. 

Our work was done for the case of weak lensing, but since it has only used the mathematical structure of the shear and the con- 
vergence, it applies also to other 2D polarization fields such as the polarization fluctuations of the Cosmic Microwave Background. 
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Appendix A: Fourier transform of the F and Gq kernels 

The F and Go kernels, as defined in O and (|4]l relate the 2- and 3-pt shear and convergence correlation functions, respectively. In 
Sect. 2 we have derived their explicit forms, see ([19} and (l22T i. According to the Fourier space relations of the shear and convergence 
statistics dSTT i and ( |32] |. one expects that Fjif- and e"^''^ are Fourier pairs, as well as -Gq/tt^ and Q^^Vi+fii+fii) ^ Here we perform the 
Fourier transforms of F/n^ and -Gq/tt^, with F and Go given in ( fT9] l and (l22b . 

The Fourier transformation of the kernel F is a 2D integral of the form 



Fio _ r 

J 



(A.l) 



The Fourier transformation of the kernel Go can be greatly simplified by performing translational shifts to the integration variables. 
It turns out that the full transformation is composed of 2D integrals similar to that in (lA.lb . 



Performing the 2D integrals in polar coordinates results in 
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Table B.l. Forms of the convolution kernel l-l as defined in ( IB. lb for different weak lensing 2-pt statistics. 
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Combining dA.lb and (IA.3b yields 

which demonstrates that the Fourier transformation of Fjip- is indeed the phase factor e'*''' in (ISTb . 
For the Go kernel we still need to take account that t^, - ^3 e'^' = -t\ - (2, so that 



4- 



_ ^1^3 (^1 + ^3) ^2^3 (^2 + ^3) _ t\lltT, I J_ _ ^1^2^3 _ 2i(/3,+/J2+A) 

~ /*/* ~ /•2 I /• /* ; f*f*f* ~ 

Y '2_ 3 ^1 2' 12 3 

Thus the Fourier transformation of -Gq/tt^ equals the phase factor e^'*'''^^^^^^^ in (l32l i. as expected. 



Appendix B: Relations between otiier correlation functions 

In Sect. I2.2l we have derived the relations between the shear and the convergence 2- and 3-pt coiTelation functions. The method 
we used is based on the relation between the 2-pt coiTelation functions of the convergence k and the deflection potential ip (fT2l i. Thus 
the method can easily be generated to derive relations between the coiTelation function of k and that of any weak lensing quantity g 
which can be expressed as derivatives of if/. We denote g = Dgi/r, and write these 2-pt relations in a general form 

(^(-^i)^'(-^2)> ^^j'^^y tv) f^ixi -X2-y). (B.l) 
Listed below are some candidates for g and the corresponding operator Dg, where a is the deflection angle defined as a = dij/. 
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Let DgDg' act on both sides of the relation between (iA(A) and {kk) (fTTI ). and use the statistical homogeneity of the k field, one 
can obtain an integration form of the kernel "H, in analogy to ([8]l. Let the same operator act on both sides of (fT2l i. one can express 

as derivatives of the convolution kernel T in (fT2b . as l-lixi - X2- y) - DgDg';F(JCi - X2- y)- Further inserting the explicit form 
of T ( fT9l ) allows one to obtain the explicit form of •?/ as a function of z - X\-X2- y- We summarize some of the 2-pt relations in 
Table. IBTI Note that the form of I-ICF) for {aa) has a minus sign, which is due to the fact that dxidx2'F{x\ - X2 - y) - -d^T(z) 
with z - xi- X2 - y- 

We write the relations between {kkk) and the 3-pt correlation functions of the g's also in a uniform convolutional form, 

{gg'g"){xi,X2)^^ j dV j d^y2 {KKK)(yi,y2)I{xi-yuX2-y2)- (B.2) 

To find the explicit form of the convolution kernel /, we first write it into an integral form, in analogy to (|6), and then split it into 
terms which can be expressed also as derivatives of the kernel f, like ( l20t . Then with the explicit form of T' one can reach the 
explicit form of /. We list the forms of the convolution kernel / for some 3-pt statistics in Table. IB. 21 

Some of these relations, e.g. those for (Ky), {jjk), and {jkk), can find their application in galaxy-galaxy(-galaxy) lensing. Some 
other relations, e.g. those for {aa), (aK), and (aaK), are potentially of interest for studies of the lensing efi'ects on the Cosmic 
Microwave Background and its cross-correlation with galaxy weak-lensing maps (.1^1:^,20001) . 
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Table B.2. Forms of the convolution kernel / as defined in ( IB.2I ) for different weak lensing 3-pt statistics. 

(XXX) integration form of / split form of / I(a, b) 
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